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Abstract. - We rigorously derive the Ornstein-Zernike asymptotics of the pair-correlation 
functions for finite-range Ising ferromagnets in any dimensions and at any temperature above 
critical. 



The celebrated heuristic argument by Ornstcin and Zcrnikc 1 implies that the asymptotic 
form of the truncated two-point density correlation function of simple fluids away from the 
critical region is given by 

Gif)^^=e~^, (1) 



where the value of the inverse correlation length £^ depends only on the density p, the inverse 
temperature (3 and the spatial dimension d. The original OZ approach hinges on the assump- 
tion that the so called direct correlation function C(-), which is de facto introduced through 
the renewal type relation 

G(f) = C(r) + p f C(r- f^G^dn, (2) 

JR d 

is of an appropriately short range. 

Because of the physical significance of both the conclusions and of the underlying heuristic 
assumptions a number of works (see e.g. |2I.''>I llojj] ) were devoted to attempts to put the the- 
ory on a rigorous basis, that is to derive JQ) directly from the microscopic picture of intermolec- 
ular interactions. Most of these works, however, were based on expansion/perturbation tech- 
niques and required technical low density or high/low temperature assumptions and, thereby 
addressed the situation when the parameters are far away from the critical region. Since the 
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OZ theory was, above all, intended to describe the phenomenon of critical scattering it would 
be of interest to devise such rigorous approach to which would rely only on qualitative 
features of noncriticality such as, for example, finite compressibility (or finite susceptibility in 
the context of ferromagnetic lattice models). 

In this Letter, we present a fully nonperturbative derivation of the direction dependent 
analog of for finite-range ferromagnetic Ising model above the critical temperature in any 
dimension. Let J = {J v } V £Z d i be a collection of nonncgative real numbers such that J v = J_„ 
and J v = if \v\ > R, where R is some finite number. The (formal) Hamiltonian is then of 
the form 

H ( a ) = ~l E J y-xa- x (r y . (3) 

Our approach pertains to the high temperature region (3 < (3 C = /3 C (J, d), which is the set of 
all (3 such that the susceptibility x/3 = ^2x£Z d ( a o a x)/3 ls fi^c- By the Simon-Lieb argument 
Xf3 < oo implies strict exponential decay of the two-point function. Equivalcntly, the series 

Xf)(t) = E e (M W,>/3, (4) 

where (■, ■) is the usual scalar product in M. d , has a nonempty domain of convergence for 
each subcritical value of the inverse temperature (3 < f3 c . Note that, by an important result 
of Aizenman, Barsky and Fernandez [7], (3 C is actually the usual critical temperature: The 
spontaneous magnetization is positive whenever (3 > (3 C . 

In the sequel we shall use to denote the domain of convergence of (0J. From a purely 
geometric point of view the direction dependent inverse correlation length £p is the support 
function of Kg. In particular, the dependence of £p(n) on the direction n S § d_1 is encoded 
in the geometry of d~K.p. 

Theorem 1 In any dimension d > 1 and for any ferromagnetic model (jSJ the asymptotic 
decay of the two-point correlation function in the high-temperature region (3 < /3 c (J,d) is 
given by 

(-oa x ) fJ c-^=e-^^(l + o(l)), (5) 
V\ x \ { ' 

where n x = x/\x\ £ S d_1 is the unit vector in the direction of x, and the function ^> p is strictly 
positive and analytic. Moreover, the inverse correlation length ^pin) is an analytic function 
of the direction n 6 S^ -1 in the sense that the boundary dK.p of is locally analytic and 
strictly convex. Furthermore, the Gaussian curvature up of cJK^ is uniformly positive, 

R/3 = min Kp(t) > 0. 
te9K,3 

Alternatively, strict convexity and analyticity of the direction dependent inverse correlation 
length £/3(-) could be formulated in terms of the geometry of the unit sphere dXJp in the £p- 
norm. 

A full proof of Theorem ^ can be found in [S]. 

In principle our approach pertains to any model in which the pair-correlation function g{x) 
admits a suitable graphical representation of the type 



A: Q^x 



(6) 
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where the sum runs over a family of paths or, more generally, pathlike objects connecting 
and x, possibly with compatibility constraints, e.g., some form of self-avoidance; let us call 
such paths admissible. The weights q( ■ ) arc supposed to be strictly positive and to possess a 
variation of the following four properties: 

• Strict exponential decay of the two-point function There exists C\ < oo such 
that, for all x £ 1 d \ {0}, 

g(x) = J2 < Cx e-«*> , (7) 

A: O^x 

where £(x) = — lim/ ; ^ 00 (fc|x|)~ 1 log <j([fcr]) is the inverse correlation length. 

• Finite energy condition: For any pair of compatible paths A and rj define the condi- 
tional weight 

q(X\i]) = q(\llr))/q(r)) 

where A II 77 denotes the concatenation of A and ry. Then there exists a universal finite 
constant C2 < 00 such that the conditional weights are controlled in terms of path sizes 
|A| as: 

q(X\v) > e~ C2 ' A| - (8) 

• Splitting property: There exists C3 < 00, such that, for all x, y £ Z d \ {0} with x 7^ y, 

E lW< C -2 E 9(A) E 9(A). (9) 

A: — >x — >y A: — >x A: x — >y 

• Exponential mixing : There exists C4 < 00 and 8 £ (0, 1) such that, for any four 
paths A, 77, 71 and 72, with A II r\ II 71 and A II r] II 72 both admissible, 

4^ll<exp{C 4 E °™}- (10) 
g(A|77ll72) ^ 

^e7iU72 

Many models enjoy a graphical representation of correlation functions of the form JHJ. In 
perturbative regimes, cluster expansions provide a generic example. Nonpcrturbative exam- 
ples include the random-cluster representation for Potts (and other) models 9 , random line 
representation for Ising |10II11II12| , or more generally random walk representation of N- vector 
models [13) . etc... However, it might not always be easy, or even possible, to establish prop- 
erties , © , © and (fLU)) for the corresponding weights. 

Before proceeding to explain how such an expansion is used in order to prove Theorem ^ 
let us first discuss the similar, but much simpler case of self-avoiding walks (SAW), which has 
been treated, though with a somewhat different approach, in [141115] . 

Self-avoiding Walks. Here we are interested in exact asymptotics of the following quantity: 

g^(x) = E , 

A: 0— >X 

where the sum runs over all finite-range SAW connecting and x, |A| is the length (i.e. the 
number of steps) of A, and (5 > 0. It is known ^B] that gj| AW is well-defined and, accordingly, 
that © holds for all (3 as soon as (3 > f3f AW (d). The remaining three properties trivially 
follow with C 2 = /?, C 3 = 1 and C* 4 = 0. 
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KVg+Xl 





KUp + xo 



Fig. 1 - The construction of the ^-skeleton \k = (xo, ■ ■ ■ , xn) of a path A = (zo, . . . , z n ). Set xo = zo; 
then define iteratively Xk+i to be the first point of A outside the set KTJp + Xk- When this procedure 
stops set xm = z n . 



We are now going to decompose a path A : — > x into elementary irreducible pieces. The 
notion of irreducibility is adjusted to the direction of the target point x through the geometry 
of the inverse correlation length £| AM/ : Let K^ AW/ be the convex set supported by £| . 
Choose a dual point t € dK^ AW such that (t, x) = £@ (x). We say that a vertex i G A is a 
regeneration point of A if the hypcrplanc through i, orthogonal to t, cuts A into two pieces. A 
path is then said to be irreducible if it does not contain any regeneration points. Since SAW- 
weights possess the factorization property e _/3 ' 7UA ' = e - ' 3 ' 7 ' e _/3 ' A ', we arrive to the following 
Ornstein-Zcrnikc type equation, or in the probabilistic jargon a renewal equation: 

g r w (x) = Y.4 AW (x-y)9r w (y)> (u) 

2/GZ d 

where the direct two-point function is defined via the summation of path weights over irre- 
ducible paths, 

c% AW (y)= T e-^l. 



■-A r (y)= E 

X-.O^y 
ftrcduciblc 

The short range nature of c^ AW finds then its precise mathematical expression in the claim 
that the direct correlation length is strictly smaller than l/^ AW (x) : 

liminf-(fclxl)- 1 logcl AW ([fcx]) > ^ AW (^) ■ (12) 

k— >oo 

It is very easy to establish l|12ll at large values of (3. The main difficulty is to give a 
nonperturbative proof. As one gets close to /3f AW , paths typically have a very complicated 
messy structure at the lattice scale, and one expects them to behave properly only at distances 
large compared to the correlation length. It is thus natural to introduce a coarsegrained 
description of the microscopic paths. To this end, we choose some big number K, and construct 
the K -skeleton \k of a path A as described in Fig.^ Notice that KXJp are balls of radius K 
in the metric £| AW , so both the scale and the geometry of the inverse correlation length ^ 
literally set up the stage for our path coarsegraining procedures. Given a target point x we 
use a dual direction t € dK^ AW ; (t,x) = £p(x), to measure the amount of backtracking by 
different skeleton steps: For any v G Z d define the surcharge cost s t {v) = £,p AW (v) — (t, v) > 0. 
Accordingly, the surcharge cost of a skeleton Xk : ^ x is defined as sum of all the surcharge 
costs of all its steps, St(Xx) = J2k 5t ( Xk ~ ^fe-i)- 

The following energy-entropy argument explains what is gained by the above change of 
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scale: By the splitting property the weight of a skeleton Xk ■ i— > x is bounded above as 
q^ AW (X K )= ]T e ~ m <e-^ A "»-^). 

A^ \k 

On the other hand, the total number of all AT-skeletons of length N is of the order (K^ 1 ) 1 ^ = 
exp{c AHog If}. Therefore, energy dominates entropy as soon as K is large. In particular, up 
to exponentially small probabilities, only a negligible fraction of Xk -steps arc backtracking 
(have surcharges of order K) with respect to x. Because of the finite energy condition (JSJ) 
such predominant forward structure of paths on the fixed finite scale K implies the massgap 
property l(T2j l. 

Repeated iterations of give: 

oo n 

^(^"^E E Il c I AW ^) e(t,!/<) 

n—1 yi+y2 + ---y n —x i—1 

which, in probabilistic terms, is tantamount to independence of different irreducible paths. 
Since c^ AW (y) e^' y ^ is a probability measure on 7L d when t G dK.^ AW , the proof of the OZ- 
formula for self-avoiding walks reduces to a local limit computation for sums of independent 
random variables with exponentially decaying tails. A very similar construction applies to all 
subcritical short range Bernoulli bond percolation models [171118] . 

Ising Model. A convenient graphical representation for the pair-correlation functions of the 
Ising model is the random-line representation |11U12|ITU| . 

gp {x) = (a x <J y ) fj = 9/3 (A) , 

A: x^>y 

where the sum runs over admissible paths connecting x and y, making jumps only between 
sites u and v such that J u -v > 0, and qp is some strictly positive weight. It can then be 
proved that properties J7J, ©, © and (fTUfl hold for this representation. We would like to 
proceed similarly to what was done for SAW. The major new difficulty arising now is that for 
the Ising model the weights do not factorize: qp(X II 7) 7^ qp (A) 9/3(7). The closest expression 
to {TJ one can write now is 

^W=E E 5/3 (A) E 9/3(7 1 A). (13) 

y£Z d A: y — >x T-O^y 

irreducible 7IIA admissible 

The admissibility constraint does not create serious problems, for example it is automatically 
satisfied across the regeneration points. However, the presence of the conditional weight 
9/3(7 I A) instead of 9/3(7) destroys the independence and, consequently, dramatically changes 
the probabilistic structure of the equation. 

At this point one is compelled to adjust the traditional point of view on the nature of 
the OZ-formula and to try to understand it not through the spectral theory of the renewal 
type equations @ or (|11|) . but rather in a broader context of local limit properties of one- 
dimensional systems of objects (irreducible paths) under appropriate mixing conditions. As 
we shall explain below such an approach leads to a reformulation of the problem in terms of the 
statistical mechanics of Ruelle's operator for shifts over countable alphabets. The OZ-formula 
ijSJ is recovered in this way from a local limit analysis based on the analytic perturbation 
theory of the corresponding Perron-Frobenius eigenvalue. 
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For the moment let us postpone the exact definition of path-irrcducibility and focus on 
the induced irreducible representation of gp{x): 

9?( x ) = E E 9,3 (Ax II... II A„). (14) 

n>l A: O^x 

A=AiU...UA„ 
Afe irreducible, fc—l,...,n 

Repeated iterations of suggest to rewrite the weights as 

n-1 

qp(\i II ... II A„) = qp(X n ) Y\_ <?/3(Afc | A^+i II • ■ ■ II A„) 

fe=i 

Introducing a dummy empty path and using it to extend finite strings of irreducible paths 
(Ai, . . . , A n ) to infinite strings A = (Ai, . . . , A ra , 0, 0, . . . ), we rewrite <|14[) in a more uniform 
and compact form: 

^0*0 = E E cxp{V(r fc A)}, (15) 

n>l \:0i-^x 

where r is the shift tA = r(Ai, A2, . . . ) = (A2, . . . ) and the potential tp is defined via e^-' = 
qp(\i IA2 II ... ) under the convention qp(\ | II II ... ) = g^(A). 

The expansion 1|15|) suggests to introduce the normalized operator L as follows: 

Lf(X)= E e^+^^/^A), (16) 

v irreducible 

where V(y) G 1> d is the displacement along the path u, and, as before, t G dUp satisfies 
(t,x) = £p(x). 

()15() falls into the framework of the classical Gaussian local limit theory once the potential 
ip happens to be Holder continuous: Given two different strings A = (Ai, A2, . . . ) and A = 
(A' l5 A 2 , . . . ) of irreducible paths define i(A, A') = min{fc > 1 : Xk 7^ A' fe }. We would like to 
find 6 G (0, 1) such that 

,,, . , . f \m-^)\ ^ ri7 s 

\ip\g = mf ml —. '- < 00. (17) 

fc i(A,A')>fc k 

Clearly, the exponential mixing property (|10|) alone is not sufficient to ensure l|17|) : Even if 
fc > 1 the sum J^«eAi ^2vex k S^ v ~ u ^ is beyond control unless one imposes further restrictions 
on the geometry of irreducible paths. This is precisely the motivation behind the following 
refined definition of irreducible paths: Given t G <9K/3 and a rcnormalization scale K, let us 
say that i\ is a break point of a path A = (ii, . . . , i n ) if ii is a regeneration point of A and 

{ii+i, . . . ,i„} Cii+ KVp +Cs(t) , 

where C$(t) = {v G Z d : (t,v) > (1 — 3){;p(v)} is a positive cone along the direction of t. A 
path A is irreducible if it does not contain break points. If Ai II • • • II A„ is the concatenation 
of such irreducible paths, then for every 1 < I < k < n, 



which, in view of (|10fl . already implies that that the potential ip in (|16|l satisfies the Holder 
condition (|17|l . It remains to check that up to exponentially negligible weights typical paths 
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7 : t— > a; contain a density of break points. However, coarsegraining procedures similar to 
those developed above for SAW, and based on the properties Q , 10 and © , show that there 
exist M < oo and v > such that 

irreducible 

This is the analogue of the mass separation result (|12fl . 

Both the local limit asymptotics of Theorem ^ an d the geometry of can be now read 
from the analytic dependence of the leading eigenvalue p(z) of L z on the perturbation z E C d 
in 

LJ(X)=L(e^»f) . 

Indeed, log p(z) is nothing but the limiting log- moment generation function of V(Xi). 
On the other hand, since is the domain of convergence of Q}, the equation of dUp in a 
neighborhood of t is given by {t + z E R d : p(z) = 1}. In other words, locally the surface 
dJip is just a level set of p through t. Furthermore, choosing perturbations z of the form 
z = p\n x + ip, where p lies in the tangent hyperplanc to &K.p at t, we recover the results of 
Paes-Leme ( Proposition 4.2) on the analyticity of the 1-particle mass shells. 

* * * 

M.C. and Y.V. gratefully acknowledge the kind hospitality of Technion where part of this 
work was done. 
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